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Abstract 



We prove the existence and uniqueness of a viscosity solution of the parabohc 
00 ■ variational inequality (PVI) with a mixed nonlinear multivalued Neumann-Dirichlet 

I boundary condition: 



"^^^^^ - Ctu {t, x) +dip{u{t, x)) 3 f{t, X, u{t, x),{Vua){t, x)), t > 0, x eV, 

^^Q^^^ + di;{u{t, x)) 3 g{t, x, u{t, x)) , t > 0, x e Bd {V) , 
n(0, x) = x gV, 

where dip and 5-0 are subdifferentials operators and Ct is a second differential 
operator given by 

d'^v{x) , / dv {x) 



C-tv{x) = \Y.{aa*),,{t,x)^^ + Y,h{t,x) 



2 ^ dxidxj ^ dxi 

i,j=i ■' 1=1 

The result is obtained by a stochastic approach. First we study the following 
backward stochastic generalized variational inequality: 

r dYt+F {t, Yt, Zt) dt+G (t, Yt) dAt G dip (Yt) dt+di; (Yt) dAt+ZtdWt , < t < T, 
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where (^t)t>o is a continuous one-dimensional increasing measurable process, and 
then we obtain a Feynman-Kag representation formula for the viscosity solution 
of the PVI problem. 

AMS Classification subjects: 35D05, 35K85, 60H10, 60H30, 47J20, 49J40. 

Keywords: Variational inequalities; Backward stochastic differential equations; 
Neumann-Dirichlet boundary conditions; Viscosity solutions; Feynman-Kag formula. 

1 Introduction 

Viscosity solutions were introduced by M.G. Crandall and P.L. Lions in and then 
developed in the classical work of M.G. Crandall, H. Ishii, P.L. Lions [2], where are 
presented several equivalent ways to formulate the notion of such type solutions. The 
framework of this theory allows for merely continuous functions to be the solutions of 
fully nonlinear equations of second order which provides a very general existence and 
uniqueness theorems. 

In 1992 E. Pardoux and S. Peng [10] introduced backward stochastic differential 
equations (BSDE) and supplied probabilistic formulas for the viscosity solutions of 
semilinear partial differential equations, both of parabolic and elliptic type in whole 
space. Elliptic equations with Dirichlet boundary condition have been treated by 
R.W.R. Darling and E. Pardoux in [3] and with a homogeneous Neumann boundary 
condition by Y. Hu in [1] . 

The parabolic (and elliptic) systems of partial differential equations (PVI without 
the sub differential operator) with nonlinear Neumann boundary conditions was the 
subject of the paper E. Pardoux and S. Zhang ^A]. The case of the systems of vari- 
ational inequalities for partial differential equations in whole space was studied by L. 
Maticiuc, E. Pardoux, A. Ra§canu and A. Zalinescu in [7j. 

The main idea for proving the existence of the viscosity solutions for PDE and 
PVI is the stochastic approach. Using a suitable BSDE, or backward stochastic vari- 
ational inequality (BSVI) for the PVI case, one can obtain a generalizations of the 
Feynman-Kag formula (i.e. stochastic representation formula of the viscosity solution 
for deterministic problems). 

The origin of our study goes from the PDE 

— -Ctu = f,t>o,xeV, 
' -^ = 0, t>o, xeBd{v), 

^ u{0,x) = h{x), X eV, 

which is a mathematical model for the evolution of a state u{t, x) G M of a diffusion 
dynamical system with sources / acting in the interior of the domain T> and g on the 
boundary of V. 
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In certain applications it is call upon to maintain the state u{t, x) in a interval 
I C M for all X & V and in a interval J C M for all x G Bd{T>). Practically these 
can be realized adding the supplementary sources dh {u {t, x)) and dl^ {u {t, x)) on the 
system. These sources produce "inward pushes" that would keep the state process 

u{t, x) in I,WxeV and u(t, x) in J, Va; G Bd (V) 

and do this in a minimal way (i.e. only when u{t, x) arrives on the boundary of I and 
respectively J). Hence dli{u (t,x)) and dij {u{t,x)) represent perfect feedback flux 
controls. 

The aim of this paper is to treat the more general case of a parabolic variational 
inequality with mixed nonlinear multivalued Neumann-Dirichlet boundary condition. 
This requires the presence of a new terms in the associated BSVI considered, namely 
an integral with respect to a continuous increasing process. 

The scalar BSDE with one-sided reflection, which provides a probabilistic represen- 
tation for the unique viscosity solution of an obstacle problem for a nonlinear parabolic 
PDE, was considered by N. El Karoui, C. Kapoudjian, E. Pardoux, S. Peng, M.C. 
Quenez in [3]. E. Pardoux and A. Ra§canu in f8j (and [9] for the generalization to the 
Hilbert spaces framework) studied the general case of BSVI and obtained probabilistic 
representation for the solutions of PVI in whole space. 

The paper is organized as follows: In Section 2 we formulate the Neumann-Dirichlet 
PVI problem; we present the main results and we prove the uniqueness theorem. For 
the existence theorem we first study in Section 3 a certain BSVI. The solution of 
this backward equation gives us, via Feynman-Kag representation formula, a viscosity 
solution for the deterministic multivalued partial differential equation as it is shown in 
Section 4. 



2 Main results 

Let V he a. open connected bounded subset of M"^ of form 

V = {x eM.'^ : £{x) <0} , Bd iV) = {x G M'^ : ^ (x) = O} , 

where i G Cl (M^), \V^ {x)\ = 1, for all x e Bd (V). 
We define outward normal derivative by 

^ - E HV. (.) . V. (.)) . for all . . (V) . 

j=l 3 3 

The aim of this paper is to study the existence and uniqueness of a viscosity solution for 
the following parabolic variational inequality (PVI) with a mixed nonlinear multivalued 
Neumann-Dirichlet boundary condition: 
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du{t, x) 
dt 

du{t, x) 
On 



— CtU (t, x) + d(f[u{t, x)) 3 f(t, X, u{t, x), (VMCr)(t, x)), 

t>0, xev, 

+ dtlj{u{t,x)) 3 g{t,x,u{t,x)), t > 0, X e Bd{V), 



(1) 



u(0, x) = h{x), X eV, 
where operator Ct is given by 



Ctv{x) = -TT[a{t,x)a*{t,x)D^v{x)] + {b{t,x),Vv{x)) 



dv{x) 
dxi 



for V eC^ [W^) . 

We will make the following assumptions: 



(I) Functions 



h : [0, oo)xW^ ^ W^, 
a : [0, oo) X M"' ^ M"'^'^, 
/ : [0, cx)) X ^ X M X M'^ ^ I 
g : [0, oo) X 5rf (D) X M ^ I 
h -.V are continuous 



(2) 



We assume that for all T > there exist a G M and L, (3,'y > (which can 
depend on T) such that Vt G [0, T] , Vx, x G M'^ : 



|6 (t, x) — b{t,x) I + ||cr (t, x) — CT (t, x) II < L ||x — x| 
and Vt G [0, T],yx eV, u e Bd (V), y,y eR,z,z e M'^: 



(3) 





{y - y) {f{t, X, y, z) - /(t, X, y, z)) < 




{^^) 


\f{t,x,y,z) - f{t,x,y,z)\ < /3\z - z 


1 


in) 


|/(t,x,y,0)| <7(1 + ||/|), 






{y - y) {git, u, y) - g{t, u, y)) <a\y- 


-y\\ 


iv) 


\9{t,u,y)\ <-i{l + \y\). 





(4) 



In fact, condition (jH-i) and (jH-if ) mean that, for all t G [0, T], x G m G Bd {V), 
z G R'^: 

y ^ ay - f {t,x,y,z) -.R 
y ^ ay - g{t,u,y) -.R^R 

are increasing functions. 



(II) With respect to functions (y9 and ip we assume 

{i) M ^ (— cxD, +00] are proper convex l.s.c. functions, 

{ii) f{y)>(f (0) = and (y) > (0) = 0, V ?/ G M, 

and there exists a positive constant M such that 



(5) 



< M, Vx G 



< M, Va; G 5rf (P) . 



(6) 



Remark 1 Condition ^ii) is generally realized by changing problem (QP in an equiv- 
alent form, as example: if {uq, Uq) G dip we can replace ip (u) by ip {u + uq) — (p (mq) — 
{uq, u); a similar transformation one can do for ip. 

We denote 

Dam {(p) = {m G M : (/? (u) < 00} , 

dip iu) = {u* : u* {v - u) + ip iu) < ip iv) , Vt; G M} , 

Dam {dip) = {m G M : ^v? (m) ^ 0} , 

(m, u*) e dip ^ u E Domdip, u* G dip (u) 

(for function ip we have the similar notations). 
In every point y G Dam (ip) we have 

dipiy)=Rn [ip'_iy),ip'4y)], 

where ip'_{y) and ip\{y) are left derivative and, respectively, right derivative at point 

y- 

(III) We introduce the compatibility assumptions : 

for all e > 0, t > 0, X G 5rf ("D), X G ^, y G M and 2 G M'^ 



{i) Vipe (y) g (t, X, y) < [Vipe (y) g {t, x,y)]^, 
(ii) Wipe iv) f (t, X, y, z) < [Vipe (y) f (t, x,y,z)]^, 



(7) 



where a+ = max{0,a} and Vipe{y), Vipeiy) are unique solutions U and l^, 
respectively, of equations 

dip{y -eU)3U and dip{y - eV) 3 V. 
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Remark 2 A) Clearly, using the monotonicity ofVtfejVipe, we see that, if 

y ■ 9it,x,y) <0 and y ■ f {t, x,y, z) < 

for allt>0,xE Bd(T>), x E V, y E M and z G , then compatibility assumptions 
^ are satisfied. 

B) If (fjip : ^ (~oo, +00] are convex indicator functions 

,ifye[a, 00) 



and 



where a < < b, 
then 



^ (y) = i[a,oo) (y) 



^ (y) = h-oo,b] (y) 



+00, ify^ [a, 00) 

, tfy e (-00, fe] 
+00, ify^ (-00,6] 



Vv^e (y) = ~iy- «) (^nd V^e (y) = ^ (z/ - 

and the compatibility assumptions become 

9 (t, X, y) > 0, for y < a, and 
f {t,x,y,z) <0, fory>b. 

We shall define now the notion of viscosity solution in the language of sub- and 
super-jets, see [2]. SR'^^'^ will denote below the set of d x d symmetric non-negative 
real matrices. 

Definition 3 Let u : [0, 00) x D ^ M a continuous function, and {t, x) G [0, 00) x V. 

We denote by x) (the parabolic superjet of u at {t,x)) the set of triples 

{p,q,X) G M X M.'^x^'^^^ which are such that for all {s,y) G [0,oo) x V in a neigh- 
bourhood of (t, x): 

y) < u{t, x) + p{s - t) + {q,y - x) + 

+ ^{^{y - x),y - +o{\s-t\ + \y-x\'^). 

Similarly is defined V'^~u{t, x) (the parabolic subjet of u at {t, x) ) as the set of triples 
{p,q,X) G M X R'^xSR'^^'^ which are such that for all {s,y) G [0, 00) x V in a neigh- 
bourhood of {t, x): 

u{s, y) > u{t, x) + p{s - t) + {q,y - x) + 

+^{^iy -x),y-x) + o{\s -t\ + \y- xp), 

where r — * o (r) is the Landau function i.e. o : [0, oo[ — > M zs a continuous function 
o (r) 

such that lim = 



We can give now the definition of a viscosity solution of the parabohc variational 
inequality ([1]). We denote first 

V (t, X, p, g, X) =^ p - ^Tr x)X) - {bit, x), q) 

- f{t,x,u{t,x),qa{t,x)). 

Definition 4 Let u : [0, oo) x D ^ M a continuous function, which satisfies u{0,x) = 
h{x), \/ xeV, 

(a) u is a viscosity suhsolution of IHj if: 

u(t,x) e Dom{ip) , \/{t,x) e {0,oo) xV, 
ult, x) e Dam {^) , V(t, x) E (0, oo) x Bd {V) , 

and, at any point (t,x) G (0, cxo) x V, for any {p,q,X) E V'^'^u{t, x) : 
V (t, X, p, q, X) + ip'_ {u{t, x)) <0 if xEV, 

min ^^V {t,x,p,q,X) + ip'^{u{t,x)), {Vi (x) ,q) - g{t,x,u{t,x)) (g) 

+^'_ {u{t, x)) } < tf xE Bd (V) . 

(b) the viscosity supersolution of ^ is defined in a similar manner as above, with 7^^'+ 
replaced byV^~ , the left derivative replaced by the right derivative, min by max, and 
the inequalities < by > . 

(c) a continuous function u : [0, oo) x V is a viscosity solution of (Qp if it is both a 
viscosity sub- and super- solution. 

We now present the main results 

Theorem 5 (Existence) Let assumptions ([^-(0) be satisfied. Then PVI (QP has a 
viscosity solution. 

For the proof of the existence we shall study a certain backward stochastic general- 
ized variational inequality (then we use a nonlinear representation Feynman-Kag type 
formula). We present this approach in the following section and after then the proof 
of Theorem [5] in Section |H 

Theorem 6 (Uniqueness) Let the assumptions of Theorem\^ be satisfied. If function 
r — > g{t, X, r) is decreasing for t > 0, x G Bd (T>) , (9) 

and there exists a continuous function m : [0, oo) [0, oo), m (0) = 0, such that 
|/(t,x,r,p) - /(t,?/,r,p)| < m(|x -?/| (1 + IpD), 

V t > 0, X,?/ G p G 

then the viscosity solution is unique. 
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Proof. It is sufficient to prove the uniqueness on a fixed arbitrary interval [0, T] . 

Also, it suffices to prove that if -u is a subsolution and f is a supersolution such that 
m(0, x) = f (0, x) = h{x)^ X eV, then u < v. 

Firstly, from definition of V, there exists a function i & [W^) such that i (x) > 
on V with Vi{x) = Vi (x) for x e Bd iV) (as example i{x)=l (x) + sup |£ {y)\). 

For A = |a;| + 1 and 5, £, c > let 

u (t, x) = e^*u {t, x) — 6i (x) — c 

V (t, x) = e^^v {t, x) +Si{x) +C + e/t. 



Denote 



/ (t, X, r, g, X) = \r- -Tr [ (aa*) (t, x)X] - x), q) 



f (t, X, e '^V, e '^^gcr (t, x) ) 



(11) 



and 



^(t,x,r) = e^^g{t,x,e ^V) 



Clearly r ^ f {t, x, r, q, X) is an increasing function for all {t, x, q, X) G [0, T] x 
j^d^gR'^x'^. Moreover, since 



M= sup (x) \ + \D£ (x) \ + \D'i (x) \ + \b{t, x)\ + \a{t, x)\ \ < oo, 

then we can choose c = c{6, M) > such that for u = u{t, x) and i = i (x) : 

/ (t, X, M, Dm, D^m) < /(t , X, m + 5£ + c, Dm + (5D£, + 5D'^£) . 

Using these properties, assumption (Q, and the fact that left and right derivative of 
(p,ijj are increasing we infer that function u satisfy in the viscosity sense 

— (t, x) + / (t, X, u {t, x) , Du (t, x) , D'^u {t, x) ) 



+ e^V- (e^^*M(t, x)) < if X G P 



f du 



+e^V-(e"^*'«(^,a;)) , ( V£ (x) , Dm (t, x) ) + 5 
-g{t, X, u{t, x)) + e^V- (e^^*M(t, x)) I < ii x e Bd (V) . 



(12) 



Analogously we see that v satisfy in the viscosity sense: 
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dv ~ 

^ (t, x) + f{t, X, V (t, x) , Dv (t, a;) , D^v {t, x) ) 



+ e^V+ (e"^*t;(t, a;)) - e/t^ > iixeV 



dv 



max < — (t, x) + / (t, X, f (t, x) , i^f (t, x) , (t, x) ) 
+e^V- (e~^*t;(t, x)) - e/t^ , ( V£ (x) , Dv (t, x) ) - 5 
(t, X, t;(t, x)) + e^V+ (e"^*w(t, x)) I > iixe Bd (V) 



(13) 



For simplicity of notation we continue to write u, v for u, v respectively. 
We assume now, to the contrary, that 

max_ {u — v)^ > 0. 

[0,T]xI> 



(14) 



Exactly as in Theorem 4.2 in [8J we have (t, x) G [0, T] x Bd{T>) where (t, x) is the 
maximum point, i.e. 

■u(t, x) — f (t, x) = maj^ {u — f )"*" > 0. 

[0,T]xC 

We put now (see also the proof of the Theorem 7.5 in Crandall, Ishii, Lions [2J) 
$„ (t, x,y)=u (t, x) -V (t, y) - p„ (t, x, y) , with (t, x, y) G [0, T] x D x P, 

where 



Pn {t,x,y) = ^ |x - + 5((t,x,M(t,x))(V£(x) ,x - y) + |x - 



2 



"^V- (e~^*M(t, x)) (V£ (x) , X - y). 



(15) 



Let it be (t„, x„, ?/„) a maximum point of $„. 

We observe that m (t, x) — f (t, x) — |x — x|^ — |t — has in {i, x) a unique maximum 
point. Then, by Proposition 3.7 in Crandall, Ishii, Lions [2], we have, as n — > oo 



u {tn, x„) ^ u(i, x), (t„, x„) i;(t, x). 

But domain V verify the uniform exterior sphere condition: 

3 ro > such that ^(x + roV£ (x) , ro) n "D = , for x G 5d (V) 

where S'(x, tq) denotes the closed ball of radius tq centered at x. 
Then 

|y — X — roVi (x) I > Tq , for X G Bd (V) , y E V, 



(16) 



or equivalent 



{V^ {x) ,y - x) < ^\y - for a; G Bd (V) , y e V. 

If we denote 

B {t, X, r, q) = (V£ (x) , q) - g{t, x, r) 
then, if x„ G Bd (V), we have, using the form of p„ given by f|T5|) and f|T7|) . that 

+g{i,x,u{i,x))Vi{x) + 4 |x„ - (x„ - x) - e^V^ (e"^*M(t, x)) V£ (x) 
> |x„ - +^(t,x,M(t,x))(V£(x) , V£(x„) ) 

(^rt) -^nj (^nj -^n) ) ~l~ 4 \Xn x| ^ V£ (Xji) , Xji X^ 

-e^V^ (e-^*M(t, x)){Vi (x) , V£ (x„) ) 
Then ([16]) imphes for x„ G 5rf (P) : 

> 



hm inf 

n^oo 



^gAt„^/_ (e~-^*"M(t„, x„)) 
Analogously if ?/„ G Bd (V) we infer 

limSUp B itn,yn,V {tn,yn) , -^yPn (^n, 3;^, - 5 

n—>oo L 

Then from (fT2!) . (fT3|l we conclude that 

P + f{tn, Xn, U (tn, X„) , D^pn {tn, ^n-, Vn) , X) + e^^"ip'^ (e"^*"M(t„, X„)) < 0, 



for {p,DxPnit 

ri) ^ri) y-n ),X) gP ' M(t„,xO 



and 



P + /(^n, 1/n, (^n, 2/n) , -DyPn {tn, Xn, Vn) , Y) 

+e^*'V'+(e-^*"t;(t„,?/„)) > 



t2 



2 

for (p, -DyPn (t„, x„, ?/„) , F) G "P ' v{tn, Vn) 

From Theorem 8.3 in Crandall, Ishii, Lions [2j (apply with k = 2, Oi = O2 
ui = u, U2 = —V, bi = p, &2 = —p ) we deduce that there exists 



(p, X, r ) G M X 



^dxd ^ ^dxd 
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such that 



and 



),x) eV '^uit 

{P, -Dypn {t ),Y) eV v{tn,yn) 



-(^+11^11) 



/ 
/ 



< 



X 0\ ^ , 1 ,2 

< A + -A, 
n 



-r 



where A = Dl (t„, ^n^yn)- From (lisp we have 



^2 = 2n2 



where |0 < C \h\ (the Landau symbol). Then fl20l) become 



7 
.0 / 



< 



— (3n + K 
where k„ — > 0. Now from f|T8|) and f|T9l) 



X 
-F 



< 3n 



/ -/ 
-/ / 



/ 
/ 



t2 



< f{tn,yn,v{tn,yn),-DyPn{t +e^*"(^V(e-^*"t;(t„,?/„)) 

/ (^rt) -^n, It (tnj -^jj) , DxPn ij^m •^ri) y-n) ; ^) ^ "V-'— (c "uitn, 

By (dH) and ([H]) there exists N >1 such that 

X„) > V{tn, yn), W n> N 

and consequently 
and 

f (^n; ^ (^n? •^n) ; Dypn (^n; "^nj 1/n) ? ^) ^ 

/ i^ni y-nj ^ (^n; Z/ra) y ^yPn (tny -^ny y-n) y ^) 

Then, by definition (fTTl) of / and assumption (ITOl) . we have 

-oo Z' ('^n; l/n? ^ ('-'ri) -^n) ; DyPn (tny -^ny yn) y F) 
' f (tny -^ny ^ (j^ny •^n) y DxPn (j^ny -^ny yn) y ^) 



— < liminf 
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But from 1^, ^ q,qe 

{Xq, q) - {Yq, q) < 3n |g - q\^ + ( \q\^ + 

Hence 

Tr [ {aa*) (t„, x„)X - {aa*) 

= ^ (Xo-(t„, x„)ej, a(t„, x„)ei) - {Ycr{tn,yn)ei, a{tn,yn)ei) 
1=1 

<3C n \Xn - Unf + ( \(T(tn,Xn)f + W {tn, yn)f ) l^n, 

and consequently 

^<0 

that is a contradiction. 
Then 

u (t, x) <v (t, x) , V (t, x) G [0, T] X V. 

■ 

3 Backward stochastic variational inequalities 

Let {Wt : t > 0} be a c?- dimensional standard Brownian motion defined on some com- 
plete probability space (f2,jF, P). We denote by {J-'t : t > 0} the natural filtration 
generated by {Wt : t > 0} and augmented by Af the set of P- null events of JF: 

jr^ = a{Wr : < r < t} V7V. 

Let r : Q —>■ [0, oo) be an a.s. J?-t-stopping time and 

• {^t : t > 0} be a continuous one-dimensional increasing progressively measurable 
stochastic process (p.m.s.p.) satisfying Aq = 0. 

We shall study the existence and uniqueness of a solution (Y, Z) of the following 
backward stochastic variational inequality (BSVI) : 

dYt + F (t, Yt, Zt) dt + G (t, Yt) dAt e dip (Yt) dt + dijj {Yt) dAt + ZtdWt, 

< t < r, (22) 
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3.1 Assumptions and results 



Let A, yU > 0. 
Let 



'l[o,^] (s) ds®dF; 
the Hilbert space of p.m.s.p. / : X [0, cx)) ^ such that 



E 



e \j[s)\as 



1/2 



< oo 



and 



the Hilbert space of p.m.s.p. / : X [0, oo) ^ such that 



n 



1/2 



< OO . 



We also introduce the notation S^''^ for the Banach space of p.m.s.p. f : Qx[0, 
l'^ such that 



oo 



e(sup e^*+'^^*|/(t) 



1/2 



< OO . 



With respect to BSVI ( l22l) we formulate the following assumptions: 



Let F : X [0, oo) x 



Sfe TUikxd 



G : X [0, oo) 



X 



satisfy that 



there exist a, /3 G M, L > and r], 7 : [0, 00) x ^ [0, 00) an p.m.s.p. such that 
for all t>0,y,y'e R'', z, z' E R'''"^ : 



and 



(«) F{-,-,y,z) is p.m.s.p., 

(ii) y — > F {uj, t, y, 2;) : ^ R^ is continuous, a.s. 

(m) {y - y', F (t, y,z)- F {t, y', z)) <a\y- y'\^ , a.s. 

{iv) \F {t,y, z) — F {t,y, z') \ < L\\z — z'\\ , a.s. 

{v) \F{t,y,z)\<r]t + L{\y\ + \\z\\), a.s. 

(0 G{-,-,y) is p.m.s.p., 

{ii) y — > G {(jj, t, y) : R'' ^ R^ is continuous, a.s. 

(m) {y - y', G (t, y) - G (t, y')) < f5\y - y'f , a.s. 

{iv) \G{t,y) I < 7t + L li/l , a.s. 



(23) 



(24) 
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Terminal date ^ is an M^-valued jFT-measurable random variable such that there 
exists X > 2a + 2L^ + 1, fi > 2(3 + 1 : 



Jo 



(25) 



oo. 



(26) 



Let it be ip, ip such that 

{%) ifjip : M.^ {—oo, +oo] are proper convex l.s.c. functions, 

(tt) ip{y)>ip{0) = 0, ^(y)>^(0)=0, 

The sub differentials are defined by 

dif (x) = {v e R'' : {v,y - x) + If (x) < if (y) , G M''} 
and similar for ip. 



The existence result for (!22|) will be obtained via Yosida approximations. Define 
for e > the convex C^-function if,, by 



fe (y) = inf 1^ |y - vf + ipiv):vE M'^j 



y - Jey 



(and similar for ip^). 
Denoting 

Jey ={I + edif)'^ (y) and Vfe (y) 
Hence y —>■ Vfe {y) is an monotone Lipschitz function and 

1 2 

(y) = ^\y^ •^^y\ + ^ (-^ei/) • 

(analog for ip^). 

• We introduce now the compatibility assumptions : 
for all e > 0,t>0, y eR'' and z E R'''"^ 

(z) (V<^,(2/),V^,(y)> >0, 
(tt) {Vfe {y) ,G{t,y))< ( V^, (y) , G (t, y) 
(m) (V^, (y) , F {t, y,z))< {Vfe {y) , F {t, y, z) )" 



(27) 
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Definition 7 (Y, Z, U, V) will be called a solution of BSVI ^ if 

(a) Y G n nt'^ n n',^, z e nit,, 
{h) u e nl'^, V e nl'^, 

(c) E / e^^^'^^^ {Ys) ds + tp (r,) dAs) < oo, 



(d) {Yt, Ut) e d^, F {du) ® dt , {Yt, Vt) e dip, F {duo) A (cj, dt) 

a.e. on Q X [0, r] 

(e) Yt+ [ Usds+ [ VsdAs = ^+ [ F{s,Y,,Zs)ds 



(28) 



tAT 



tAT 



tAT 



+ G{s, Ys) dAs - / ZsdWs, for all t > a.s. 

JtAr JtAr 



In all that follows, C denotes a constant, which may depend only on /i, a, j3 and L, 
which may vary from line to line. 

Proposition 8 Let assumptions and If {Y, Z, U, V) and {Y, Z, U, V) 

are corresponding solutions to C, and C, which satisfy ^25\) . then 



Jo 



\Ys - {ds + dAs) + \\Zs- ZsW^ds 



+E sup e^^+^''^'\Yt -Yt\^ <CE 

0<t<T 



(29) 



Proof. From Ito's formula we have 
gA(tA.)+M.A.|y^^^_y^^^|2^ T e''+^^^\Ys - Yf {\ds + fldA^) 



tAT 



+2 / 6^^+^^^= - Us - Us)ds + 2 / 6^^+^^^= {Ys - Ys, Vs - Vs)dAs 

It At J tAT 



e^'+^'^^WZs- ZsW'ds 



tAT 



^gAr+M.|^_||2 + 2 / e'''+^''^{Ys-Ys,F{s,Ys,Zs)-F{s,Ys,Zs))ds 

JtAT 

+2 / e^'+''^^{Ys-Ys,G{s,Ys)-G{s,Ys))dAs 

JtAT 

-2 / e^^+^^^(r, - Ys, {Zs - Zs)dWs) 



tAT 



Since 



{Ys - Ys, Us - Us)ds > 0, {Ys - Ys, K - Vs)dAs > 0, 
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2(n - F (s, F„ Z,) - F{s, Zs)) < (2a + + - + -| - 
and 

2(n-n,G(s,F,)-GGs,n)) < (2/? + i)|n-y;|2 

then (using also the Burkholder-Davis-Gundy's inequahty), inequahty fl29|) follows. 
The main result of this section is given by 



Theorem 9 Let assumptions ^23\)-(21) he satisfied. Then there exists a unique solu- 
tion (r, Z, U, V) for 



3.2 BSVI - proof of the existence 

Consider the approximating equation 

y/ + r Vv^, (r/) ds + r v^, (r/) dA = ^+ I f (s, r/, zt) ds 



tAT 



iAr 



tAT 



+ G{s, y/) dAs - / zidWs,yt > 0, p 

JtAr JtAr 



(30) 



a.s. 



Since Vv?£, V'i/'e : M'^ ^ M'^ are Lipschitz functions then, by a standard argument 
(Banach fixed point theorem when y F {t, y, z) and y — » G (t, y) are uniformly 
Lipschitz functions and Lipschitz approximations when y ^ ay — F {t, y, z) and y 
j3y — G(t,y) are continuous monotone functions) (see also pjj), equation (I5UI) has a 
unique solution 

(F^ z^) G (5^^ n Hi'^ n H^'^) X n'^i. 



Proposition 10 Let assumptions ^^23^-^20^ he satisfied. Then 



E 



sup 

0<t< 



r Jo 



Xs + fiAs 



+ 1 e 

'0 



(31) 



< C M (r) 



Proof. Ito's formula for e'*'*+^"^* l^tf yields 



iAr 



+2 / e^"+^^= 



(r/, v^, (r/) )\ds + (F/, v^, (y/) )ixdA, 



+2 / e^^+'^^^(F/,F(s,y/,Zf))c/s + 2 / e^^+^^= (F/, G (s, F/) 



tAT 
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f T 



-2 / 6^^+^^^= {Y^.ZldW,). 

JthT 

But from Schwartz's inequality and assumptions (l23ll -(|25ll we obtain 

2(n, F (s, Z,) ) < 2a + 2L \Y,\ \\Z,\\ +2\Y,\ \F {s, 0, 0) | 

< (2a + 2L2 + 1) |r/ + i ||Z,f + |F (s, 0, 0) |' 

and 

2(y;, G (s, r„ z,) > < 2/3 |r/ + 2 |y;| (s, o) | < (2/3 + 1) |y;|' + 1^(5, 0) 

Hence, using also that (y, Vv?^ (y)) , (y, V^/^j {y)) > 0, 

itAr 

+ / 6^'+'^^= (\F{s, 0, 0) I'rfs + |G (s, 0) l^dAs) 



Hl\T 



that clearly yields (for A > 2a + 2L2 + 1 and > 2/? + 1): 

E r e'^^"'"' [ [ds + dA,) + ||Z,^f ds\ <C M (r) 
Since, by Burkholder-Davis-Gundy's inequality, 

\ 1/2 



E sup 

0<f<r 



< ^E sup e^*+^^' + CE r e'^'"' \\Zl\\' ds 

4 0<t<r Jo 



then it follows 

E sup e^*+^^' < C M (r) 

0<t<T 

The proof is complete. 
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Proposition 11 Let assumptions ^2&^-^^^ he satisfied. Then there exists a positive 



constant C such that for any stopping time 6 G [0, t] : 

(a) e^^+'^^^(^|Vy5e (^/) \^ds + |VV^, (F/) \^dA,^ <C M{t) , 

ih) E re"^+^^^((^(j,(r/))ds + ^(j,(F/))rfA,) < C M{t) , 

(32) 

(c) Ee^^+^'^« - J, (r/) f + |r/ - J, (y/) < ^ c m (r) , 

(d) Ee^^+^'^« (^(J, (Yi) )+^{Je (F,-) )) < C M (r) . 



Proof. Essential for the proof is the stochastic sub differential inequality introduced 
by Pardoux and Ra§canu in [8], 1998. We will use this inequality for our purpose. First 
we write the subdifferential inequality 



for s = tj+i A T,r = ti At, where t = to < ti < t2 < ... < t A t and tj+i — ti = —, then 

n 

summing up over and passing to the limit as n ^ oo, we deduce 

JtAT ^ 

+ [ ^e{Y!)die'^^^^^) 



We have the similar inequalities for function ipir 

If we summing and we use equation fl30l) . we infer that for all t > : 

gA(M.)+M... (^^ (y^e^j + (y^. J ) + r e^^+^^» I V^. (F/) 

Jt^T 

JtAr JtAr 

+ f e^^+'^^=(v^, (r/) , v^. (r/) > {ds + dA,) 

< e^^+^^^ (v., (0 + i^e (0 ) + e^^+^^^(V^, (F/) , F (s, F/, Zf) >ds 



'tAr 
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+ 



JthT 



The result follows by combining this with ( 13T1) . assumptions ( 1271) and the following 
inequalities 

1 I ,2 1 I ' ,2 

( Vy., {y) , F > < i I Vy., (y) |' + 3(r^2 + ^2 |^|2 ^ ^2 1 1^| |2 ^ 

(V^, (y) , G y) > < i I V^. + 2(7,2 + \yf ) 
( (y) , F (s, y,z))< (Wife (y) , F (s, y, z) 

(Vyp, {y) , G (s, > < (VV^, (y) , G (s, 

<^ivv^.(z/)r +2(7.^+^2 11/1') ■ 

Proposition 12 Lei assumptions ^2B^) -^^2l\j he satisfied. Then 
E r e'*^"''(|y/-y/|'(ds + rf^) + ||Zf-Zf||'t/s) 

+ E sup 6^*+'^^* |y/ - y/|' < G (e + (5) M (r) 

0<i<T 

Proof. By Ito's formula 
gA(tAr)+M*A. I y^e ^ _ ^ 1 2 ^ /" gA.+M. | y^e _ 1 2 ^^^^ ^ ^^^^^ 

JtAr 

Jt^T 

+2 / 6^^+^^^= (F/ - y/, (r/) - VV^5(F/)> rfA, = 
= 2 r e^^+^^^ (r/ - y/, F (., F/, - F(., y/, )> 

JthT 

+2 r 6^^+^^^= (r/ - y/, G (s, y/) - g{s, y/)> rfA, 



(33) 



'Mr 
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- [ 6^^+^^^= \\Zl - Zf\\^ ds-2 [ 6^'+"^^ (F/ - F/, (Z| - Z^^)dWs) 

JthT JtAr 

eover, 



'iAr 

We have moreover 



<(2a + 2L2)|r/-F/|^ + l||Zf-Z 



But from definition of and the monotonicity of operator dip we have 

= (V^, (F/) - VMYs'), Y! - r/> - e| v^, (r/) f - 6\v^s{y!)\^ 
+ {e + 5){w^e (y/),v^5(r/)) 

then 

(Vy., (y/) - V¥.,(y/), F/ - F/> > - (5 + 5) (Vyp, (F/) , V^s{Y^)) 
and in the same manner 

(vv'e (n^) - VMY^). y: - r/> > - (5 + 5) (v^. (r/) , vv^^ir/)) 

and consequently 



JtAT 

JtAT 2 

<2{e + 6) r e'^^^^'\{Vve{Y!),VMY!))ds 
Now, from (l^a) 



tAT 



2{e + S)E e 



tAT 



+ VMYf))dA, <C {e + 6) M {r) 



(34) 



and clearly by standard calculus inequality (1331) follows. ■ 

We give now the proof of Theorem M 
Proof. Uniqueness is a consequence of Proposition [81 The existence of solution 
(F, Z, U, V) is obtained as limit of (F/, Z^ Vipe , V^-e (F/)) 
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From Proposition [12] we have 



limY' = Y in S^''' nUl''' nH^ , 

e\0 



hmZ^ = Z in Ti. 

e\0 



kxd' 



Also, from (l32]-a) and (l32]-c') we have 



hmj,(r") = r in nl''', hmj,(r^) = r in 



\,fi 



£\0 



£\0 



hmEe^^+^^* I Js{Yf) - = 0, hmEe^^+''^'' 1 7^(1^) - = 



£\0 



£\0 



for any stopping time 9, < 6 < t 

Using Fatou's Lemma, from fl32] -b). fl32] -d) and the fact that is l.s.c. we obtained 
(IMc). 

Denoting = V<^e(Y'), V = Vi^eiX') , from (I32l-a) it follows: 



E 



6^'+"^' {\U'\'^ds + \ V'\^dA,) 



<C M{t) 



Hence there exists U G and V G TY^''^ such that for a subsequence 

jje-r, _^ weakly in Hilbert space Ul'^ 
ye„ weakly in Hilbert space Ti^''^ 



and then 



E 



< lim inf E 

n— >oo 



{\Ufds+\VfdA,) 



Je 



<C M2 {9, t) 



Passing now to lim in (130|) we obtain (!28] -e). 

Let u G nl'^, V G Ht'". Since V<^e(>"/) e dip{Je{Yf)) and V^/',(r/) G 
Vt > 0, then as signed measures on x [0, r] 

gA.+MA,^f/£^ - (F/) ) P (duj) 0ds + e^'+>'^^ip{Je (F/) ) P (dcu) ds 

< e^"+^^= P (dio) ® ds 

and 

^gAs+M.^(j^ (r/) ) P (rfcj) A (to, ds) 

< e^'+f'^^'^ivs) P (du) ® A {u, ds) . 
Passing to lim inf in these last two inequalities we obtain fl28] -d). The proof is complete. 
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4 PVI - proof of the existence theorem 

It follows from result in [6] that for each (t, x) G M+ x V there exists a unique pair of 
continuous JF*-p.m.s.p. (X*'^, AI'^)s>q, with values in D x M_)_, solution of the reflected 
stochastic differential equation 

rsVt psVt ps\Jt 

s I — > Al'^ is increasing (35) 

psVt 

<. Jt 

where 

J^l = a {Wr - Wt : t < r < s} y Af. 
Since P is a bounded set, then 

sup I X*'^' I < M (36) 

s>0 

and with similar calculus as in [11] we have for all n,T,p > there exists a positive 
constant C such that V t, t G [0, T] , x,x' : 

E sup \X'f - Xf <C{\x- x'\P +\t- (37) 
se[o,T] 

and 

Ele^^^T] < oo. (38) 

Let T > be arbitrary fixed. Under assumptions (l2])-(I7j), it follows from Theorem 
M with r replaced by T that for each {t, x) G [0, T] x V there exists a unique solution 
(F*^, Z*^, f/*^', l/*^') of p.m.s.p. 

Z*^ G H^'^, t/*^ G T^i'^, V*^ G Ki'^ 

with y/'^ = r/'"^, Z*'^ = 0, f/*'^ = 0, 1/*'^ = 0, for all s G [0, t] 
solution of BSDE: 

J s J s 

+ ri[*,T](r) f{T,X'/X'''.Zl^^)dT 

J s 

+ /" \t,T] (r) (7 (r, X*'", F,*'") dA^'" - /" for all s G [0, T] a.s. 
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such that (F/'^, f/*'^) G d^, P {duj) ® dt, (y/-^, V;*-^) G SV', P {duj) ® A {uo, dt) , a.e. on 
Vt X [t, T] . 

We observe that function /, g depends by uj only via function X*'^. 

Proposition 13 Under assumptions we have 

E sup 6^'+^"^^ I F/'^' I ^ < C (T) (39) 

se[o,T] 



and 



E sup e^^+^^^lF^'^ - < E e^"+'^^^|/i(X^^') - h{Xtr''')\^ 

sG[0,T] L 

+ re^^+^^^|l[,T](r)/(r,Xf ) - l[,,,^](r)/(r, )|'tir (40) 

^0 



T 



+ / e^^+^^'■|l[^,T](r)^?(r,X*^y;*-•) - ^,^T]{r)g{r,X';-' X^fdA'^^ 



Proof. Inequahty ( l39l) follows from Theorem [9] using also ( l36l) . ( l38l) . Inequality f j40ll 
follows from (l29l) in Proposition [HI ■ 

We define 

«(^,x) = F/^ (t,a;) G [0,T] xP (41) 

which is a determinist quantity since Y^^ is JF* = A/'-measurable. 
From Markov property we have 

n(.,Xf)=rf (42) 
Corollary 14 Under assumptions {^-(0), function u satisfies: 
(a) M(t, x) G L'om ((^) , V(t, x) G [0, T] x ^, 

(6) M(t, x) G Dom (V^) , V(t, x) G [0, T] x Bd (V) , (43) 
(c) M G C([0,T] x^). 

Proof. Using (l28l -c) we have (f(u{t,x)) = E,Lp(Y^^) < +oo and similarly for Hence 
(jl3]-a.b) follows. Let (tn,Xn) — > (t,x). Then 

|M(t„,x„) - M(t,x)|' = E|F/;"" - < 2Esup,g[o,T] in*""" - 

+2E|r/^ - F/^P 

Using (jlU]) . fl5Bl) . fl?r|l and (1551) we obtain M(t„,x„) u(t,x) as (t„,x„) — > (t,x). ■ 
We present now the proof of Theorem [S] (existence of the viscosity solutions). 
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Proof. It suffices to show the existence of the solution of PVI ([T]) on an arbitrary fixed 
interval [0,T]. Setting 

u {t, x) = u (T — t, x) 
then the existence for problem ([T]) it is equivalent with existence for 



+ Am {t, x) + f{t, X, u{t, x), {Vua){t, x)) G d(p{u{t, x)) , 

t e {o,T) , X ev, 
du{t,x) ^ g(j-^^^^(^^^^^'^ g dtp{u{t,x)), t e (0,T) , X e Bd{V) 



(44) 



where 



and 



dn 

u(T,x) = h{x), X eV 



f {t, X, u,z) = f {T - t, X, u,z), g {t, x,u) = g{T - t, x, u) 
a (t, x) = a (T — t,x) , b (t, x) = b(T — t, x) 



d^v ix) ~ , dv (x) 



XI 



2 dxidxi ^ dxi 

We denote also 



. . l\JJ-.r. 



V^(t,x,p,g,X) _p_ iTr((aa*)(t,x)X) - {h{t,x),q) 
- f{t,x,u{t,x),qa{t,x)). 

In the sequel, for simplicity we keep notations b, a, u, f, g, C, V instead of b, a, it, /, g, C, V 
and we shall prove that function u defined by (|4T1) is a viscosity solution of parabolic 
variational inequality (jS]). We show only that m is a viscosity subsolution of (14^ (the 
supersolution case is similar). 

Let (t, x) e [0, T]xV and (p, q, X) G P^'+u{t, x). 

1. The proof for the G P is similar of that from [8]. 

2. Let X G Bd{V). Suppose, contrary to our claim, that 

min (t, X, p, g, X) + ip'_ (u{t, x)) , 

( V£ (x) ,q) -g {t, X, u{t, x)) + {u{t, x)) I > 
and we will find a contradiction. 
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It follows by continuity of /, u, b, a, i, left continuity and monotonicity of </?'_ 
and that there exists e > 0, S > such that for all \s — t\ < 5, \y — x\ < S, 



-ip + s)- -Tr y) {X + si)) - y),q+{X + si) {y - x)) 

-f{s, y, u{s, y),{q + {X + si) {y - x) )a{s, y)) + (fL{{u{s, y)) > 0, 



(45) 



if x e -D 



and 



{V^ {y),q+{X + si) (y-x))- g{s, y, u{s, y)) + y)) > 0, 

ifxeBd (V) 

Now since {p, q, X) e V'^'~^u{t, x) there exists < S' < S such that 

u(s,y) < u{s,y), 

for all s e [0,T], s t, y & V, y ^ X such that |s — 1| < S', \y — x\ < S', where 
u{s, y) = u{t, x) + {p + s) {s - t) + {q,y - x) + ^{{X + si) {y-x), y-x) 



(46) 



Let 

We note that 



u inf {s>t: \X'f -x\> S'} 
(^M^ ) = {Y^'\ ) ,t<s<{t + S')Ai. 



solves the BSDE 

F/'^ = u (i/, ■) + / (/(r, Z^'-) - Ut''^)dr 



+ 



J s 



t,x 
r 



(y/'^, C/*'^) e P (du;) di, (y/'^, T/*'^) e dip, P (du;) A {uj, dt) , 

a.e. on O X [t, T] . 

Moreover, it follows from Ito's formula that 

(n*'^ ) = {u{s, X'f), iyua) {s, )) ^t<s<t + 5' 

satisfies 

'du{r, Xl'^ 



Y^''^ = u{u,Xln 



i: 



dt 



+ Cru{r,Xl''') 



dr 



Zt^^'dWr 



+ / {V,u{T,X'/),Vl{Xl^''))dAY 
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Let (y;*'^, z*'^) = (i;*'^ - y/-^, zl'^ - zl'^). 

We have 



du{r, Xl'^) 



dt 



t,x 



dr- / ZTdW, 



+ 



Let 



4 = X*.-) + /(s, X*'-, y/'-, i*.-) 

Since — < C {Zf^ — Zl'^\, there exists a bounded d— dimensional p.m.s.p. 
{C; 0<s<iy} such that4 - Ps ^ (C, ^^''^) 



Now 



+ 



+ 



dt 



dr 



It is easily to see that, for the process 

1 



r* = exp 



|C.|'cir+ / {Cr.dWr) 



we have, from Ito's formula. 



and so 



dt 



t,x\ 



rfs + r*(z;'" + y;'"Cs,cii^s) 
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Then 



-E 
-E 



dt 



t,x 



dr 



(47) 



We first note that (Yi, f/t) G (F, "K^) G dip imphes that 

X'f))ds < U'fds, tpi {u{s, < V^^'^'dA 

Moreover, the choice of 6' and u imphes that 

uiiy,Xin<u{^,X'f) 

From (US]) and (USD it fohows that 

-{p + e)-ps + X*'")) > 0, if X G P 

and 



t,x 



duis.X'f) 
dn 



All these inequalities and equation (H7|) imply that F/''^ > and equivalent 

which is a contradiction with the definition of u. Hence we have 
min|v"(t,x,p,g,X) + {u{t, x)) , {Vi (x) , q) - g{t, x,u{t, x)) + ip'_{u{t, x))^ < 



This proves that u is a viscosity subsolution of (IH|) . Symmetric arguments show that 
u is also a supersolution; hence u is a viscosity solution of PVl 



Remark 15 Ifb,a,f and g do not depend on t then we have a directly a representation 
formula for the viscosity solution u of PVI 



0,x:t 



u {t, x) = Yq 

where Z^'^'*, ?7°'^'*, V;°'^'*)o<s<t solution of BSVI 



+ f ^(X°'",i;°'"'*)ciA°'"- /" Z°'"'*ciPV„ for all s e [0,T] a.s. 

J s J s 
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an 



d (Xf'^ /lO'^)o<.<t solves SDE 



Jo Jo 



Jo 



0,x\j aO,x 
r 1 



A 



0,x 



A^'^ is increasing 



1 J /\0,x 

^{X'^'='eBd{V)}"'^r ■ 



Corollary 16 We have 

u{t,x) e Dom{d(f) , V(t, x) G [0, T] X P 

Proof. Let {t, x) be fixed. We have two cases: 

1) Dom {dip) = Dom {ip) , and so, from (l43l-a). u{t^ x) G Dom {dip). 

2) Dom {dip) 7^ Dom {ip). Let b G Dom ip \ Dom {dip). 

Then b = sup{Domip) or b = inf Dom ip. If 5 = sup{Domip) and u{t,x) = b, then 
(0,0,0) G V'^'+u{t,x) since 

u{s,y) < u{t, x) + o{^\s — t\ + \y — x\'^) 

and from ([8]) it follows ip'__{b) = ip'_[u{t,x)) < oo and consequently b G Dom {dip); a 
contradiction which shows that u{t,x) < b. Similarly for b = m{{Domip). ■ 
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